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Quantum gravity "foam", among its various generic Lorentz non-invariant effects, 
■ would cause neutrino mixing. It is shown here that, if the foam is manifested as 

o : 

, a nonrenormalizable effect at scale M, the oscillation length generically decreases 



X 



with energy E as (E/M) 2 . Neutrino observatories and long-baseline experiments 



■ should have therefore already observed foam-induced oscillations, even if M is as high 

, as the Planck energy scale. The null results, which can be further strengthened by 



better analysis of current data and future experiments, can be taken as experimental 
evidence that Lorentz invariance is fully preserved at the Planck scale, as is the case 



> 
On 
O 

■ in critical string theory. 

O 

The noisy vacuum of a quantum theory of gravity could a priori be imagined to have a 
^ variety of effects on the wavefunction of a particle traveling through it. Although unitarity 

ft, ■ would probably constrain the effects of the vacuum on the wavefunction's amplitude, one 

could imagine that its phase might be shifted by the local effects of quantum black holes 



and the like which quickly pop in and out of the vacuum. Also, the type (flavor) of the 
particle might be affected: for example, a virtual or quantum black hole could as Hawking 
has suggested |1|], erase information by swallowing a neutrino of one type and spit out one 
of a different type as it evaporates. Thus, possible effects that could be considered include 
Lorentz invariance violating (LIV) indices of refraction [Q, ffl , flavor oscillations of neutral 
particles^], [7|, ffl, and corrugations of wavefronts, i.e. Rayleigh scattering M. Rayleigh 
scattering, which would imply momentum non-conservation, appears to be ruled out by the 
observation that 10 MeV neutrinos propagated without Rayleigh scattering from supernova 
1987A H (assuming that their scattering cross-section off cells of Planck foam scales as E 2 
). We conclude that if such "foam" effects exist, they must preserve translation invariance 
on the Planck distance scale and on larger ones; thus any non-vanishing effects must be 
cooperative and coherent over space-like distances in this range of scales; it may cause wave 
dispersion but not opalescence or translucence. 
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There has been considerable work on dispersive, LIV vacua. Colladay and Kostolecky 



have exhaustively classified renormalizable LIV effects [JTO], [Tl[ . Coleman and Glashow |g, [7| 
discuss neutrino oscillation from renormalizable LIV terms, which yield an oscillation length 
that scales as 1/E, and their results imply that any such effect that is of order unity at 
the Planck scale is already ruled out by many orders of magnitude. Amelino-Camelia et al. 
considered a photon velocity that has a linear term proportional to E/M and showed that 
it is marginally consistent for astrophysical gamma ray pulses if M is of order 10 16 GeV. 

The leading quantum gravity theory - critical string theory - predicts that in flat (empty) 
space all interactions are Lorentz invariant (LI) down to, and including, the Planck scale. 
LI in flat space is crucial to the internal mathematical consistency of string theory in that 
it guarantees that its symmetries are valid at the quantum level (see, for example, 
However, the Universe, possessing a preferred frame of reference, is clearly not LI. String 
theory allows such spontaneous breaking of LI, by allowing the possibility that the curvature 
and additional moduli fields have time- and space-dependent expectation values, but it does 
not allow LIV Planck scale "foam". 

In this paper we consider a particular class of LIV effects - nonrenormalizable effects 
of the vacuum, which could be induced by a high energy/short distance physical cutoff - 
and show that they give rise to a neutrino dispersion that is linear in E/M. This is the 
dispersion relation considered by Amelino-Camelia et al. M for photons, but, like Coleman 
and Glashow|B|, |7J, we consider neutrino oscillations that it could induce. In contrast to 
the latter, we consider nonrenormalizable corrections, which depend on a higher power of 
E/M. If the difference between the propagation velocities of different neutrino flavors is 
proportional to E/M then the oscillation length is proportional to (E/M)~ 2 . (Amusingly, 
for M ~ Mpi anc k it would give an oscillation length of order an astronomical unit for solar 
neutrino energies and atmospheric length scales for atmospheric neutrino energies @, but see 
below.) We then apply the argument || that high energy experiments of neutrino oscillations 
with E~ 2 mixing lengths should be able to "detect" nonrenormalizable cutoffs as high as 
the Planck scale. We argue that such experiments can, by not detecting generic effects of 
such cutoffs, provide experimental support for critical string theory versus many proposed 
alternative theories, as the latter imply LIV physical cutoffs, while critical string theory 
provides a LI cutoff. 

It is clear that our arguments are generic, and restrict or rule out only those models of 
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quantum foam in which neutrino flavor is unprotected by a symmetry and Lorentz invariance 



is broken by nonrenormalizable interaction terms, as pointed out in |H| |TJJ. It remains to 
be seen whether any specific models of quantum gravity effects can evade our conclusions, 



as suggested in |T3|, |I4f . 

To understand the appearance of nonrenormalizable LIV terms we model any type of short 
distance interaction which may be induced by quantum gravity foam, black holes swallowing 
particles and spitting them out, ultraviolet LIV cutoff physics, etc. First, we assume that 
such interaction is strong only at some characteristic energy scale M (e.g. Mpi anc k) so it can 
be modeled by the following interaction term in the effective Lagrangian, 

£ int = Mj d 4 x>^ A p (x)f* B [M(x - x')} r B (x'). (1) 

Here / determines the interaction strength, and is therefore assumed to be small when its 
argument is larger than unity, A, B = 1,2 label different neutrino eignestates of the total 
Hamiltonian, p, a — 1, ...,4 are fermionic indices. To avoid Rayleigh scattering as discussed 
previously we have assumed that / is a function only of the distance four-vector x — x'. We 
may now expand ip%(x') = ij)%{x) + d^ip^{x){x' — xY + ^d^d u ipB(x)(x' — xY(x' — x) v + 
where the dots stand for higher derivative terms. This expansion may be used to convert the 
interaction Lagrangian ([I]) into a series of local terms. Each additional derivative comes with 
an additional power of 1/M, demonstrating that the terms induced by (Q) indeed depend on 
powers of E/M. The lowest dimension nonrenormalizable terms are dimension five operators 
and therefore depend on E 2 /M. They are given by the second moments of jf B - It may 
well be that the renormalizable terms which could have been induced by ([!]) are absent. For 
example, if /4. B depends only on the absolute value \x — x'\, no dimension four operators 
are induced. If, in addition, neutrino masses are protected, say, by chiral symmetry, and 
fp B (\x—x'\) factors into constant matrix M^ B and a universal kernel f(x — x'), then neither 
are dimension three operators induced. 

Consider the diagram in Fig. 1. The effective interaction vertex connecting the incoming 
and outgoing states includes, by assumption, quantum gravity effects. Because any flavor 
gravitates, the effective vertex does not generically respect the global flavor symmetry (e.g. 
a quantum black hole swallowing one type of neutrino and evaporating into another type) so 
in its presence the eigenstates of the weak interactions \a), \/3) are mixtures of the eigenstates 
of total Hamiltonian \A), \B), \a) = cos9\A) + sm9\B), and \(3) = -sm.9\A) + cos9\B), 
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FIG. 1: Effective interaction vertex induced by quantum gravity effects. 

and therefore the interaction terms obtained from (|1|) induce mixing between flavor states. 
We will argue shortly that the induced mixing angles are generically large sin 2 (2$) ~ 1. (In 
|13| it is claimed that in certain "kinematical models" flavor mixing does not arise). 

The possibility of neutrino oscillations arises, in general, when terms are added to the 
neutrino sector of the Standard Model Lagrangian, such that the eigenstates of the total 
Hamiltonian are not coincident with, but rather linear combinations of the different neutrino 
flavors. In the simplified case of two flavors, and a single oscillation inducing interaction, 
a neutrino of flavor a can oscillate into a neutrino of flavor (3 after having traveled for a 
distance X, with a probability that is usually expressed in terms of a mixing angle and an 
oscillation length L 

-ttX 



P a ^f3 = sin 2 (20) sin 2 



(2) 



The oscillation length L = E _ Eg , is defined in terms of the difference of the states' energy 
Ea — Eb and does not depend (for this case) on the mixing angle 0. The mixing angle 
determines only the amplitude of the oscillations. In the case of more than two neutrino 
flavors, eq. (0) becomes more complicated, and includes more than one oscillatory term. 
Since eq. (|2[) contains the essential physical ingredients of the oscillation phenomenon, and 
is much simpler than the corresponding expression in the general case, we will concentrate 
on the two flavor case. 

If oscillations are induced by more than one interaction, L should be replaced by the total 
oscillation length L tot , which is given by a combination of the individual oscillation lengths 
L n that each individual interaction would have induced on its own. In the two flavor case 
(0 n are the mixing angles of each single effect) one obtains f?j] 

E^ 1 cos20 n 

Kl = ~ cos2fl tot ' (3) 
from which it is clear that in general the total oscillation length is dominated by the shortest 
one. This feature remains valid also for the case of more than two flavors, though eq. (0) is 
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replaced by a more complicated expression. 

Thus, if we are interested in a case in which one oscillation inducing term (say Lj) 
dominates, we may replace L by Lj instead of using L = L tot in eq. and use the 
factorization property of eq. (0) which implies that the oscillation length does not depend on 
the mixing angles. 

The most studied oscillation mechanism is due to the presence of neutrino mass. For this 
case, LI implies that the oscillation length for ultrarelativistic neutrinos is given by L ~ j^z, 
Am 2 being the neutrino squared mass difference and here c = 1 = h. In fact, every LI term 
in the effective action will give the same oscillation length energy dependence as a mass 
term. On the other hand, energy dependence different from L oc E is a smoking-gun signal 
of violation of LI. Several examples of LIV effects have appeared in the literature: L oc E° 
|, |, L oc E- 1 0, |TJ and also L oc E~ 2 @, |16| and L oc E' 3 0. 



Recent SuperKamiokande (SK) results indicate that a fi-r mass mixing mechanism is 
the best candidate for explaining the observed n~e anomaly in the atmospheric neutrino 



flux|T8[. The analysis of energy dependence of SK data excludes LIV effects as the primary 
source for the observed fi-e anomaly. Furthermore, by analyzing models in which LIV terms 
are present in addition to the mass mixing, it is possible to put stringent upper limits on 
their strength [|T^]. We show below that existing experimental data are already sufficient 
to provide important clues about the nature of neutrino interaction at the Planck scale, 
and that with better analysis and with more data LIV terms could be ruled out to high 
accuracy. Our claim is based on a preliminary analysis of the L oc E~ 2 case, which are 
induced only by nonrenormalizable terms in the neutrino effective Lagrangian. Dimensional 
analysis shows that renormalizable terms in the neutrino effective action induce oscillation 
lengths proportional to E~ l , or to a non-negative power of E. In particular, all possible 
renormalizable and rotation invariant mixing terms in the neutrino effective action induce 
oscillation lengths proportional either to E°, to E, or to E~ x . This can be checked by a 
case-by-case analysis of the dispersion relation induced by all the operators that are bilinear 
in the neutrino field, and have mass dimension not exceeding four. A complete set of such 
operators is ip^ip, tp^^ip, ipd t ip, xpj 5 d t ilj, ip^dtip, ip7 5 7°dtip, ip^ 0l diip, ipY 3 9 k '^ijk^ ■ (We 
have not included in this list operators that can be obtained as linear combinations of them 
and LI terms). It follows that only terms that are both nonrenormalizable and LIV can be 
the origin of L oc E~ 2 . 
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, and the total 



We now wish to show that generic mixing angles that are induced are large. Assuming 
that gravity is "flavor blind" the interaction Hamiltonian induced by the terms in eq. ([!]), 

I g{E) g{E) 

expressed in flavor basis is a "democratic" matrix H int = 

\g(E) g{E) 

Hamiltonian (including possibly mass mixing) in the flavor basis could be expressed as 

( /in + g(E) + g(E) \ 
H a p = , where h\ 2 = h 2 \- In the limit that the energy difference 

\h 21 + g(E) h 22 + g(E) J 

between the Hamiltonian eigenstates AE = \E% — E 2 \ is dominated by g{E), the Hamiltonian 

. (A , ( 1 ^ 

eigenvectors are given by -j= I , -j= I I , and therefore mixing is maximal. In case 

|<7 (-£?)! > \h\ 2 \, \hu—h 22 \ we therefore expect a large mixing angle. Since the SK (and possibly 
also solar neutrino) data appear to imply a large mixing angle due to LI effects, at least 
among some pairs of neutrino types, we infer that h\ 2 > \hn — h 22 \, so that |g(-E)| > \h\ 2 \ 
implies ^(-E 1 )! > \h\\ — h 22 \ and is by itself sufficient to infer a large mixing angle. 

As a simple explicit example of the impact of a local term resulting from the interaction 
Lagrangian ([[]), consider adding a dimension five LIV interaction term to the standard 
kinetic term of massless neutrinos, 

cAB _ 

C = ib AB ^ A $ tp B - —*pAdtda ^ B , (4) 

with S AB = diag(SA, Sb)- To obtain S AB from the second moments of f A a B we have used 
the 7-matrices as a basis to 4 by 4 matrices. In terms of the "democratic" interaction 
Hamiltonian which we have discussed previously one finds g(E) = ^jj(Sa — Sb)- 

As we have argued, an interaction term as in eq. (^) is generic to every nonrenormalizable 
theory that is also LIV. We will explicitly show that the resulting oscillation length satisfies 
L oc E~ 2 . Other terms such as jjipdida ^, jjipd t diy?p, -^^dtdt^^ip, ^dida 5 ^ ^, 
■ l j%l)dtdi r ) br )' l ip also produce L oc E~ 2 . Recall that in this case the oscillation length depends 
only on energy difference between the oscillating states, and does not depend on the mixing 
angle between flavor eigenstates which determines only the amplitude of oscillations. 

The equation of motion for the physical eigenstates is (p M 7 M — SA fj E 7°) ^a,b = ; by 
multiplying this equation by p M 7 M — Sa '^ [ E 7° and using the anticommutation relations 
of gamma matrices, we obtain E 2 (l — ^f— ) — p 2 = 0, and therefore, to first order in 
S a ,bE/M E ~ |p| (l + %f^), from which it follows that 

2tt M 



Sa-SbE*' (5) 
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independently of the mixing angle between flavor eigenstates. 

A LIV term could also arise in an otherwise LI theory, such as critical string theory, via 
spontaneous breaking of Lorentz symmetry. For example, a LIV contribution of the form 
could be generated, when a field (which could be the dilaton, or any other of the moduli 
fields of string theory) is slowly rolling in a background of a non-vanishing gravitational 
potential. But in this case, the breaking would be proportional to the time derivative of 
the field ($/M) whose magnitude is severely constrained. If $ is non- vanishing then $ has 
kinetic energy. The requirement that the kinetic energy density of $ is not larger than 
universe closure density provides an incredibly stringent bound $/Mpi anc k < 1CT 61 , which 
makes spontaneous breaking effects completely undetectable. This conclusion extends to all 
the terms that can produce an L oc E~ 2 behavior since all of them must contain at least 
one time (or space) derivative of a field, and therefore are subject to the same (or similar) 
phenomenological constraints. 

Because all nonrenormalizable theories are expected to generate higher derivative terms 
at the cutoff scale, and because such terms are not necessarily LI unless the physical cutoff 
mechanism itself is intrinsically LI, as in critical string theory, a detection (or exclusion) 
of L oc E~ 2 carries with it information about the high energy/short-distance properties 
of the theory. In contrast, the detection of an energy dependence L oc E^ 1 produced by 
renormalizable terms (though in any case extremely interesting) would not carry with it 
any such information. If LIV at the cutoff scale is explicit instead of spontaneous, typical 
mixing terms are expected to be of order unity (and detectable, as we show below), simply 
because there is no symmetry that protects them. It follows that the detection of L oc E~ 2 
effect would be a definite signal that nonrenormalizable and LIV terms are present in the 
neutrino effective action. We conclude that the detection of L oc E~ 2 dependence in neutrino 
oscillations would constitute strong evidence against critical string theory, in which the 
ultraviolet cutoff is realized in a LI way. Conversely, an experimental exclusion of such 
an effect should be considered as vindication of critical string theory as compared to other 
Planck scale models that violate Lorentz symmetry. 

If one evaluates eq.(^), with E ~ 1 GeV, M ~ M P i anck , and Sa — Sb = 0(1) as expected 
in explicit LIV where g(M) ~ M, one obtains an oscillation length of about 10 km. Typical 
high end of the energy range of solar neutrino observatories is about 10 MeV which results 
in an oscillation length of about 10 -3 AU. These estimates suggests that if the mixing 
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is close to maximal mixing sin 2 (2#) = 1, then previous neutrino oscillation experiments 
CHORUS and NOMAD |2TJ, and the currently operating experiments SKJTJ, K2KJ22 



and SN0[|3|], are already able to detect the proposed effect if Sa,b is larger than about 
1CT 2 . In this context, the existing data must be interpreted as a null result, Sa,b less than 
about 10~ 2 , or very small mixing sin 2 {26) <C 1 or both, accordingly to equation (0), as 
the analysis of the energy dependence of the oscillation length of SK data shows that the 
observed oscillation is due to mass mixing L cx E, rather than to LIV terms Planned 



experiments MINOS and CNGS[gJ| will be able to strengthen and verify these findings. 
The situation is summarized in table I which assumes maximal mixing. As can be seen 
from the table, we can already conclude without further analysis that the actual upper 
bound on a is about 10~ 3 for maximal mixing, and that the forthcoming experiments are 
potentially capable of improving this limit by at least one order of magnitude. Translating 
this bound on a into a bound on M, we obtain M > 10 3 Mpi anc k. It is worth noticing that 
the experimental sensitivity for other possible sources of nonrenormalizable LIV, such as 
the detection of time delays of photons from distant astrophysical sources, turns out to be 



significantly lower M ~ Mpi anc k (see Table I in [25]). A more detailed analysis could enhance 
the sensitivity of neutrino oscillation experiments compared to the estimates in the table, 
as shown in Fig. 2, and determine exclusion regions in a, sin 2 {26) plane. 

Furthermore, all these experiments observe (or plan to observe) a broad spectrum of 
incoming neutrino energies, from about 10 MeV to few hundreds GeV or more, so in fact 
each experiment sets a stronger constraint on nonrenormalizable LIV terms as the high 
energy end of its range is exploited. In the case of SuperKamiokande the flight-distance also 
varies - from a few tens of Kilometers to about 10 4 Km. Thus, better analysis of available 
data, and additional data from planned experiments will allow the significant strengthening 
of the bound on Sa,b, perhaps down to 1CU 8 . Figure 2 demonstrates that the limits reported 
in table I can be improved by the analysis of the high-energy part of the neutrino spectra. 
This is in contrast to the mass mixing case for which L is increasing with energy, where the 
best bounds are obtained by the lowest energy tail of neutrino spectra. Up-going muons 
in the SuperKamiokande experiments, reaching energies of 10 3 GeV , can provide the best 
bound on a at about 10~ 8 , while a value of 10~ 6 can potentially be reached by MINOS 
and CNGS. A more detailed analysis is required to determine the highest energy at which 
enough data can be accumulated. 
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EXP. 


STATUS 


(E) (GeV) 


L (Km) 


X (Km) 


a 


CHORUS 


closed 1997 


26 


io- 2 


0.85 


io- 2 


NOMAD 


closed 1999 


24 


10~ 2 


0.94 


10~ 2 


SK 


operating 


1.3 


10 


10-10 4 


1-10" 3 


K2K 


operating 


1.3 


10 


250 


10~ 2 


SNO 


operating 


0.008 


10 5 


10 8 


10~ 3 


MINOS 


starting 2003 


15 


0.1 


730 


10~ 4 


CNGS 


starting 2005 


17 


0.1 


732 


10~ 4 



TABLE I: Shown for each experiment are its operation status, mean value of observed neutrino 
energy, oscillation length according to (0) (with Sa — Sb = 1), typical neutrino flight distance X, 
and the ratio a = L/X. Parameter a can be thought of as the value of Sa — Sb for which L = X 
for each experiment, that is, the lowest value of Sa — Sb to which each experiment is sensitive, 
assuming maximal mixing. 

We conclude that existing bounds are at most barely compatible with the existence 
of a LIV ultraviolet cutoff at the Planck scale. Strengthening of these bounds by better 
analysis and additional data from future experiments would be an imminent vindication of 
critical string theory, and a strike against models that allow explicit breaking of the Lorentz 
symmetry at the cutoff scale. 
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